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Orbital degrees of freedom of a Cooper pair play an important role in the unconventional 
superconductivity. To elucidate the orbital effect in the Kondo problem, we investigated a single 
magnetic impurity coupled to Cooper pairs with a p x + ip y (d x 2_ y 2 + id xy ) symmetry using the 
numerical renormalization group method. It is found that the ground state is always a spin 
doublet. The analytical solution for the strong coupling limit explicitly shows that the orbital 
dynamics of the Cooper pair generates the spin 1/2 of the ground state. 
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Unconventional superconductivity is characterized by 
an angular momentum (p, d, or /) of its Cooper pair. It 
is very sensitive to impurities and surface boundaries, dif- 
fering from the standard BCS (s-wave) superconductiv- 
ity. As a result, it displays various phenomena associated 
with breaking of the Cooper pairs. In the last decade, 
much attention has been paid to the effect 
magnetic impurities and static boundaries Sn 
At present, the attention shifts to the study of magnetic 
impurities in the unconventional superconductivity. Re- 
cent experimental studies provided evidence of induced 
local moments with-spin S — 1/2 in Zn-doped d-wave 
superconductors.B&H A possibility of the KondOpeffect 
was indicated in impurity-doped YBa2Cu307_a.lLirE3 At 
the present stage, the search for the Kondo behavior in 
unconventional superconductors is still in progress. 

A local spin in a superconducting state causes the fol- 
lowing significant effects: (1) breaking time-reversal in- 
variance and (2) dynamical coupling with quasiparticles 
due to the Kondo effect. For the first-effect, the local 
spin can be treated in a classical wayJiJ On the other 
hand, many-body correlations play an important role in 
causing the second effect. In the case of s-wave super- 
conductors, we have a competition between the Kondo 
effect and the superconducting energy gap A. Even if the 
gap exists in the density of quasiparticle states, the local 
spin forms a spin singlet with the quasiparticles under 
the condition Ik > A, where Ik is the Kondo temper- 
ature.E-sliS A similar competition is found in a gapless 
d x 2_ y 2-wa,ve superconductor.!!!! In this case, the gapless 
density of states stabilizes the Kondo singlet. 

In this letter, we study a local spin coupled to a gapped 
superconducting state with a p x + ip y (or d x 2_ y 2 + id xy ) 
symmetry. The d- vector for the p x + ip y -~w&ve is given 
by dk = (k x + ik v )z, where z is a unit vector in the z 
axis. Recently, the property of the p x + ip y -w&ve su- 
perconductivity has been investigated extensively, which 
was motivated by the dispay^ynfthe spin triplet su- 
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perconductor S^RuO^ 



is analogous to 



the property of the A phase of the superfluid 3 He. In 
recent theoretical studies of the p x + ip y superconduc- 
tivity,|-upque physics concepts suctuas vortex charging 
effect,E3Ef spontaneous currcutEp£3 and spontaneous 
Hall effect have been reported.c3'E3) The effect of mag- 
netic impurities should also reveal rich properties pecu- 
liar to the orbital degrees of freedom of the Cooper pair. 

The Cooper pair of the unconventional superconduc- 
tivity is constructed by conduction electrons with various 
angular momenta. In the p x + ip y -w&ve state, only two 
angular momenta are coupled to each other. It enables 
us to apply the numerical renormalization group (NRG) 
met hod. £9* This is the first NRG study of such a non-s- 
wave superconducting state. Our result shows that the 
coupled conduction electron system displays a new type 
of Kondo effect. We find that the ground state is a spin 
doublet in all the parameter region of 7k/ A. The spin 
singlet ground state cannot be realized even in the strong 
coupling region (Ik /A ^> 1), since the paired electrons 
generate the spin 1/2 of the ground state. It is differ- 
ent from the s-wave case where the interchange of spin 
doublet and singlet ground states occursES We discuss 
the difference between the s-wave and p x + ip^-wave, and 
show how the spin doublet ground state is realized for 
the p x + zpy-wave in the strong coupling limit. 

In this letter, we study a single magnetic impurity at 
the center in two-dimensional superconducting systems. 
Since we treat short-range scattering here, the impurity 
couples with only the electrons having no angular mo- 
mentum (1 = 0). For the p x + ip^-wave, the total angular 
momentum of the Cooper pair is equal to one. The order 
parameter is expressed by Ae l ^ k , where <j>k is the angle 
of the Fermi vector measured from the k x axis. Since the 
angular momentum is a good quantum number for the 
Px + ipy-w&ve, the I = and I = 1 orbitals are coupled, 
and they are decoupled from the other orbitals. This is 
the reason why the Kondo problem in the p x + ip y -w&ve 
state can be treated within the two angular momentum 
spaces. Therefore, we can apply the NRG method to the 
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Kondo problem for the p x +rip v -w&ve as we do to the 



two-channel Kondo problemHii' Let us begin with the 
following Hamiltonian H — H, 



Pm+iPy 



kcr 1=0,1 



kcr 

H imp — - 2J (-JS ■ Vera' + V6 atCr ')al Qa ak'0a'- (1) 

kk' crcr' 

Here, i? Px+ i Py is the Hamiltonian for the p x + ip y -~w&ve 
pairing state, and Hi mp represents exchange interaction 
and potential scattering of conduction electrons due to 
the impurity, ej, is the kinetic energy of the conduction 
electron with wave number k, and A is the p x + ip y -w&ve 
order parameter. The subscripts Z(=0,1) and c(=T>4) 
in akicr represent the angular momentum and spin of 
the conduction electron, respectively. S and cr are 
the S — 1/2 spin operator and the Pauli matrix for 
the magnetic impurity and conduction electrons, respec- 
tively. J(< 0) and V are the antiferromagnetic and non- 
magnetic couplings, respectively. Note that the model 
given by eq. (1) is not a simple two-band Kondo model. 
Only the I = conduction electrons are coupled directly 
to the local spin. In the normal state, the I = 1 electrons 
are completely decoupled from the local spin. However, 
in the superconducting state, they are coupled to the 
I = conduction electrons via the p x + ip y -w&ve order 
parameter. We have to take the I = 1 electrons into ac- 
count as well for the p x + ip y -w&ve. We can elucidate 
the orbital effect of the p x + ip y -wtxve and compare it 
with the s-wave, since both superconductors are full-gap 
systems and have the same density of states. We can 

+ id xy -~wewe if 



obtain a similar Hamiltonian for a d x 2_ 
we replace iA by —aA and I — 1 by I — 2. 

In accordance with the Wilson's NRG procedure,Ef 
the bare conduction band is discretized logarithmically 
and the Hamiltonian (1) can be transformed intp_the 
following hopping type with staggered potentialalj'E£P 
for both p x 
states: 



ipy-wave and d x 2_ y 2 + id xy -wave pairing 



H N+l = K l ' 2 H 
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The subscript r = ± represents the two channels con- 
structed by the I = and I — 1 orbitals. Note that 
the channels are not independent due to the channel-flip 
(Ett') t crms - This means that only one of the paired 
electrons (7 = 0) couples with the local spin as in eq. 
(1), since the operator of the conduction electron at the 
impurity is written by 

E a k,l=0,cr = CAT=0,t=+,ct + CN =0yT= -^. (3) 
k 

Here, cm to is the operator of the NRG fcrmion quasipar- 
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Fig. 1. Bound state (one-hole excitation) energy (E pot /A)* mea- 
sured from the singlet ground state energy for each V. The pa- 
rameters are fixed at J = and A = 0.01. The circles are the 
NRG results. The solid line represents [1 + (qV) 2 ] -1 / 2 , where 
a = 0.929 is used to fit the NRG results. The V dependence of 
(Spot/A)* is independent of A. 



ticle in the A^-th shell. In eq. (2), £jv is given by 

e N = [i_ A -( 7V + 1 )][l-A-( 2Ar+1 )]- 1 / 2 [l-A-( 2Ar + 3 )]- 1 / 2 , 

(4) 

where A is the discretization parameter. The values with 
tilde are normalized by (1 + A _1 )/2. For simplicity, we 
used the same values for both the superconducting cutoff 
energy jand the band width, which does not alter the 
results.EJ) Throughout this letter, we keep the lowest- 
lying 500 states at each renormalization step and take 
A = 3. 

First, we discuss the case of J = to demonstrate the 
reliability of our NRG results. Figure 1 shows a finite 
V case in which an excited state (bound state) appears 
below the gap. The appearance of the bound state is 
due to the pair breaking effect of the potential scattering 
(non-magnetic impurity effect). The ratio of the renor- 
malized bound state energy and the renormalized energy 
gap converges as the renormalization step A^ increases. 
The convergence value is represented by (£^ po t/A)*. In 
Fig. 1, (_Ep t/A)* decreases with V as expected. The 
NRG result is in good agreement with the analytic solu- 
tion given by a function of VEIy 



E, 



A 



pol 



1 + {aVf 



(5) 



Let us consider the J ^ case. Since we have con- 
firmed that the potential term V does not change the 
ground state, we restrict ourselves to the V = case. 
When J — > and A ^ 0, the impurity spin is decoupled 
from the quasiparticles and there is no bound state be- 
low the superconducting energy gap as expected. In this 
case, the NRG Hamiltonian has two independent t = ± 
channels which have the same form as the s-wave ex- 
cept for the sign of the order parameter depending on r. 
The NRG energy level structure for the p x + ip y -w&ve 
has particle-hole symmetry. In the opposite limit (J / 
and A = 0) , the NRG energy level structure at the fixed 
point is given by the admixture of a strong coupling type 
and a free spin type. The former is for the / = orbital 
and the latter is for the 1 = 1. The NRG energy spectra 
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10" 3 10" 2 10" 1 10° 10 l 10' 
T K /A 



Table I. Ground and first excited states. (A) and (B) are the re- 
sults for T K /A = 2.76 X 10~ 2 and 1.38, respectively. The energy 
E is measured from the ground state (see Fig. 2). Q represents 
the total quasiparticle number measured from the half-filling. S 
is the size of the total spin, (p) and (s) are the results for the 
Px + ipy-wave and s-wave, respectively. For the p x + ip H -wave, 
the ground state is a spin doublet and the first excited state is 
a particle-hole doublet with no spin. For the s-wave, the results 
at the even number of renormalization step are given. 
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Fig. 2. Tk/A dependence of the bound state energy (B mag /A)* 
generated by the magnetic impurity. (_E mag /A)* is measured 
from the lowest-lying doublet for each Tk/A. (a) p x + ip v - 
wave case. The ground state is a spin doublet and the first 
excited state is a particle-hole doublet with no spin (see Table 
I). The results for fixed Tk are expressed by a circle and a star 
for T K = 9.20 X 10" 5 and 1.76 X 10~ 3 , respectively. Triangle-up 
and triangle-down are the results for fixed A = 0.001 and 0.005, 
respectively, (b) s- wave case. In the positive (i? mag /A)* region, 
the ground state is a spin doublet and the first excited state is a 
spin singlet. In the negative (£ mag /A)* region, the ground and 
excited states are interchanged. The level crossover takes place 
at around T K /A = 0.3.LJ 



can be described by filling quasiparticles in one-particle 
energy levels. In the A = case, an energy level lies just 
on the Fermi energy. Since this level can be occupied 
by at most two particles, the ground state has fourfold 
degeneracy. We note that the local spin is quenched 
completely for A = 0. 

When J turns on for a finite A, a bound state appears 
below the energy gap due to the magnetic impurity. Fig- 
ure 2(a) shows the Tk/A dependence of the bound state 
energy (T mag / A)* for the p x +ip y -wave. Here, the Kondo 
temperature is defined by 

T K = VPTcxp(-l/| J|). (6) 

First, we look into the small Tk/A region. As shown in 
Table 1(A), the ground state is a spin doublet. The first 
excited state is a particle-hole doublet with no spin. In 
the Kondo effect, it is favorable for the magnetic impu- 
rity to form a Kondo singlet. However, the magnetic im- 
purity has to break the Cooper pair to couple with one of 
the paired electrons. Since this costs considerable energy, 
the Kondo singlet cannot be a ground state in the small 
Tk/A region. As Tk/A increases, (E mag / A)* decreases 
monotonically and approaches the ground state energy 
asymptotically. Thus, the ground state approaches the 
A = result smoothly as Tk/A increases to infinity, 
and it becomes fourfold degenerate. This result shows 
that the spin singlet ground state is not realized in all 
the Tk/A region, which is completely different from the 



s-wave result shown in Fig. 2(b). On the other hand, 
we find that (-E mag /A)* for the p x + ip y -w&ve is scaled 
by Tk/A as for the s-wave. While the interchange of 
the ground state does not occur in the p x + ip y -w&ve 
case, almost the same Tk/A dependence of (E mag /A)* 
is obtained in both cases. For the s-wave, the energy 
gap suppresses the Kondo effect in the (i? mag /A)* > 
region, while the Kondo singlet is stabilized against 
the energy gap in the (T mag /A)* < region. ■— This 
crossover occurs at approximately Tk/A = 0.3Ji-f In 
analogy to the s-wave, there is also a crossover at around 
(-E'mag/A)* = 0.5 for the p x + ip y -w&ve. Thus, the Kondo 
effect overcomes the energy gap in the (E mag /A)* < 0.5 
region, although the ground state is still a spin doublet. 
This implies that the competition between the Kondo 
effect and the energy gap is characterized by such Tk/A 
dependence of the bound state energy for all types of su- 
perconductivity. We mention here that the NRG energy 
level structure for the d x 2_ y 2 + id xy -~w&ve is the same as 
that for the p x + ip y -w&ve. Their difference appears only 
in the wave functions. 

Let us discuss the strong coupling limit (|J| — > oo) 
case where the impurity part Hq is truncated from the 
conduction electron part. In this limit, the ground state 
property is given by Hq. Converting from the channel 
(r) representation to the angular momentum (I) repre- 
sentation 

CN=0,T=±,a — -J={fl=QM ± /i=l,o-)) (7) 

we can express Hq = iT-d + ^A, where 

Hs-d = — • CTaa'foafoa', 

oo' 

H^ = -AY J {flh* + Aoh°)- (8) 
o 

Note that the order parameter corresponds to the hop- 
ping parameter between the I = and I = 1 local orbital 
sites, since the p x + ip y -w&ve Cooper pair is formed by 
the I — and / = 1 particles. When \J\ — > oo, the I = 
particle couples with the local spin strongly to form a 
spin singlet \s) = fU |) - | Here, \a) (a =U) 
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is the wavefunction of the local spin. The spin doublet 
(fl\s)) and spin singlet (jl^fl^s), |s)) are degenerate 
since the hopping between I = and I = 1 is forbidden. 
At a finite | J|, A in eq. (8) lifts the degeneracy. The 
wavefunctions for the spin doublet (ipS=i/2) an d spin 
singlet (ips=o) are given by the linear combination of 
the following terms: 



+ C4 [It-JM) + C 5 [fl (/ T T | |) + f ol \ T) 

^=o-c 6 | s )+ C7 (/ 1 t T ||)-/ 1 t i |t)). (9) 

Here, Cj is a coefficient. The perturbation study shows 
that '0s=i/2 is the ground state, whose energy is lower 
than that of ip s =o by -||J|(A/J) 4 . When |J| is very 
large, coefficient c± of the first term in ips=i/2 is the 
largest. In the strong coupling limit, the local spin is 
almost quenched by the I — particle forming the Kondo 
singlet \s). However, the I = 1 particle is connected 
weakly to the quenched local spin \s). This does not 
mean that the I = 1 electrons destroy the Kondo singlet. 
They gain the superconducting condensation energy and 
then generate the spin 1/2 of the ground state in the 
strong coupling limit. As |J| decreases, coefficients C2 
and C3 of the second and third terms, respectively, in 
ips=i/2 increase. This indicates that the weight of the 
local spin becomes larger in the spin of the ground state. 
Thus, the ground state always has a spin 1/2 for the 
Px+iPy-wave. The Ik/ A dependence of the bound state 
energy shown in Fig. 2(a) reflects the continuous change 
of the wavefunction of the spin doublet ground state. 

Contrary to the p x + ip y -wave, the wavefunctions for 
the s-wave are given by tps=o = I s ) an d </? 3=1/2 — 
/g|/gjcr) in the strong coupling limit. In this case, the 
hopping between the I = and I = 1 orbital sites is 
replaced by the potential at the I — site: 



-A^/ot+ZoVoI- 1 )- 



(10) 



This is because the s-wave Cooper pair is formed by only 
the I — particles. The ground state is characterized by 
<Ps=o (fs=i/2) m a large (small) Tk/A region. For the 
s-wave, the I 7^ conduction electrons are completely 
decoupled from the local spin. 

In conclusion, we have studied the Kondo effect unique 
to the p x + ipy-w&ve superconductor, where the orbital 
dynamics of the Cooper pair produces the spin of the 
ground state. For large Tk, the / = electrons couple 
with a local spin strongly to form a Kondo singlet, while 
the other electrons with I = 1 are coupled to the Kondo 
singlet weakly via the superconducting order parameter. 
The NRG study has shown that the ground state is a spin 
doublet in all the Tk/A region, which is different from 
the results for the s-wave and d a; 2_j / 2-wave cases. It is 
interesting to search for a new type of Kondo effect due 
to the orbital effect of the Cooper pair in unconventional 
superconductors . 

We are grateful to H. Kusunose, R. Shiina, and K. 
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